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Abstract 

Quantum trajectory theory, developed largely in the quantum optics commu- 
nity to describe open quantum systems subjected to continuous monitoring, 
has applications in many areas of quantum physics. In this paper I present a 
simple model, using two-level quantum systems (q-bits), to illustrate the es- 
sential physics of quantum trajectories and how different monitoring schemes 
correspond to different "unravelings" of a mixed state master equation. I also 
comment briefly on the relationship of the theory to the Consistent Histories 
formalism and to spontaneous collapse models. 

I. INTRODUCTION 

Over the last ten years the theory of quantum trajectories has been developed by a wide 
variety of authord^Hl for a variety of purposes, including the ability to model continuously 
monitored open systems§'i@. improved numerical calculation!'!, and insight into the problem 
of quantum measurementol3. However, outside the small communities of quantum optical 
theory and quantum foundations, the theory remains little known and poorly understood. 

Indeed, the association with quantum foundations has convinced many observers that 
quantum trajectory theory is different from standard quantum mechanics, and therefore to 
be regarded with deep suspicion!. While it is true that stochastic Schrodinger and master 
equations of the type treated in quantum trajectories are sometimes postulated in alternative 
quantum theoriestTtJ, these same types of equations arise quite naturally in describing 
quantum systems interacting with environments (open systems) which are subjected to 
monitoring by measuring devices. In these systems, the stochastic equations arise as effective 
evolution equations, and are in no sense anything other than standard quantum mechanics 
(except, perhaps, in the trivial sense of approaching the limit of continuous measurement). 

The fact that this is not widely appreciated, even in fields which might usefully employ 
quantum trajectories, is a great pity. I attribute it to confusion about exactly what the 
theory is, and how these equations arise. 

To combat this confusion I present in this paper a simple model for quantum trajectories, 
using only two- level quantum systems {quantum bits or q-bits). Because of the recent interest 
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in quantum information a great deal of attention has been paid to the possible interactions 
and measurements of two-level systems; I will draw on this to build my model and show 
how effective stochastic evolution equations can arise due to a system interacting with a 
continuously monitored environment. 

This paper is intended as a tutorial, giving the physical basis of quantum trajectories 
without the complications arising from realistic systems and environments. As an aid in 
understanding or teaching the material, I have included a number of exercises throughout 
this paper. For the most part these involve working out derivations which arise in developing 
the model, but which are not crucial for following the main argument. The solutions to the 
exercises are included in the appendix. 

A. The model 

For simplicity, I will consider the simplest possible quantum system: a single two-level 
atom or q-bit. The environment will also consist of q-bits. While this is far easier to analyze 
than the more complicated systems used in most actual quantum trajectory simulations, it 
can illustrate almost all of the phenomena which occur in more realistic cases. Moreover, 
the recent interest in quantum computation and quantum information has spurred a great 
deal of work on such systems of q-bits, which we can draw on. 

I will assume that the environment bits interact with the system bit successively for 
finite, disjoint intervals. I will further assume that these bits are all initially in the same 
state, have no correlations with each other, and do not interact amongst themselves. (See 
figure 1.) Therefore all that we need to consider is a succession of interactions between two 
q-bits. After each environment bit interacts with the system, we can measure it in any way 
we like, and use the information obtained to update our knowledge of the state. 

This model of system and environment may seem excessively abstract, but in fact a num- 
ber of experimental systems can be approximated in this way. For instance, in experiments 
with ions in traps, residual molecules of gas can occasionally pass close to the ion, perturbing 
its internal state. In cavity QED experiments the system might be an electromagnetic mode 
inside a high-Q cavity. If there is never more than a single photon present, the mode is 
reasonably approximated as a two-level system. This mode can be probed by sending Ryd- 
berg atoms through the cavity one at a time in a superposition of two neighboring electronic 
states, and measuring the atoms' electronic states upon their emergence. In this case, the 
atoms are serving as both an environment for the cavity mode and a measurement probe. 
Similarly, in an optical microcavity the external electromagnetic field can be thought of as 
a series of incoming waves which are reflected by the cavity walls; photons leaking out from 
the cavity can move these modes from the vacuum to the first excited state. (See figure 2.) 

B. Plan of this paper 

In section 2, I describe systems of q-bits, and the possible two-bit interactions they can 
undergo. In section 3, I examine possible kinds of measurements, including positive operator 
valued and weak measurements. 
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In section 4, I present the model of the system and environment. In section 5 I show 
how it is possible to give an effective evolution for the system alone by tracing out the 
environment degrees of freedom, producing a kind of master equation. 

In section 6, I consider the effects on the system state of a measurement on the environ- 
ment, and deduce stochastic evolution equations for the system state, conditioned on the 
random outcome of the environment measurements. 

In sections 7-8, I examine different interactions and measurements, and how they lead 
to different effective evolution equations for pure states — a kind of stochastic Schrodinger 
equation. I show how under some circumstances the system appears to evolve by large 
discontinuous jumps, while for other models the state undergoes a slow diffusion in Hilbert 
space. I further show that by averaging over the different measurement outcomes we recover 
the master equation evolution of section 5. 

In section 9, I consider measurement schemes which give only partial information about 
the system. These schemes lead to stochastic master equations rather than Schrodinger 
equations. 

Finally, in section 10 I relax the assumption of measurement, and discuss how quan- 
tum trajectories can be used to give insight into the nature of measurement itself, using 
the decoherent (Consistent) histories formulation of quantum mechanics. I also briefly con- 
sider alternative quantum theories which make use of stochastic equations. In section 11 I 
summarize the paper and draw conclusions. I present the solution to the exercises in the 
appendix. 



The simplest possible quantum mechanical system is a two- level atom or q-bit, which 
has a two-dimensional Hilbert space Ti.2- There are many physical embodiments of such a 
system: the spin of a spin-1/2 particle, the polarization states of a photon, two hyperfme 
states of a trapped atom or ion, two neighboring levels of a Rydberg atom, the presence or 
absence of a photon in a microcavity, etc. All of the above have been proposed in various 
schemes for quantum information and quantum computation, and many have been used in 
actual experiments. (For a good general source on quantum computation and information, 
see Nielsen and Chuang0 and references therein.) 

By convention, we choose a particular basis and label its basis states |0) and |1), which 
we define to be the eigenstates of the Pauli spin matrix a z with eigenvalues —1 and +1, 
respectively. We similarly define the other Pauli operators a x ,a y ; linear combinations of 
these, together with the identity 1, are sufficient to produce any operator on a single q-bit. 

The most general pure state of a q-bit is 



A global phase may be assigned arbitrarily, so all physically distinct pure states of a single 
q-bit form a two-parameter space. A useful parametrization is in terms of two angular 
variables 9 and 0, 



II. TWO-LEVEL SYSTEMS AND THEIR INTERACTIONS 
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where < 6 < tt and < <p < 2ti. These two parameters define a point on the Block sphere. 
The north and south poles of the sphere represent the eigenstates of o~ z , and the eigenstates 
of a x and a y lie on the equator. Orthogonal states always lie opposite each other on the 
sphere. 

If we allow states to be mixed, we represent a q-bit by a density matrix p; the most 
general density matrix can be written 

p = p\il>){il>\ + (\-p)W)®\, (3) 

where and \ip) are two orthogonal pure states, = 0. The mixed states of a q-bit 

form a three parameter family: 

P = (^cos 2 (#/2) + ^sin 2 (#/2)) |0)(0| + 

(^sin 2 (£/2) + i^cos 2 (#/2^ |1)(1| + 

rcos(0/2)sin(0/2)(e i *|O)(l|+e-^|l)(O|) , (4) 

where 9, <p are the same angular parameters as before and < r < 1. The limit r = 1 is the 
set of pure states, parametrized as in (fj), while r = is the completely mixed state p = 1/2. 
Thus we can think of the Bloch sphere as having pure states on its surface and mixed states 
in its interior; and the distance r from the center is a measure of the state's purity. It is 
simply related to the parameter p in (H): p = (1 + r)/2, 1 — p = (1 — r)/2. 
For two q-bits, the Hilbert space Ti 2 <8> TC2 has a tensor-product basis 

AB ' 
AB ' 
AB ' 

AB i (5) 

similarly, for N q-bits we can define a basis {\iN-1iN-2 • • • ^0)}; h = 0, 1. A useful labeling of 
these 2 N basis vectors is by the integers < j < 2 N whose binary expressions are iN-i • • • ^o : 

AT-l 

\j) = \i N -i---i ), J = E ^ • ( 6 ) 

fc=0 

All states evolve according to the Schrodinger equation with some Hamiltonian H(t), 

Over a finite time this is equivalent to applying a unitary operator U to the state 

U = T :exp(-^ [/ dtH(t)} : . (8) 



l°>A* 


5|0> B 


= |oo) 


|0>A« 


5|1> B 


= |01> 


|1>A« 




= |10) 


|1>A« 




^ 111) 



to 



where T : : indicates that the integral should be taken in a time-ordered sense, with early to 
late times being composed from right to left. For the models I consider in this paper I will 



4 



treat all time evolution at the level of unitary transformations rather than explicitly solving 
the Schrodinger equation, so time can be treated as a discrete variable 

\il n ) = UnUn-l---Wo) ■ (9) 

For a mixed state p, Schrodinger time evolution is equivalent to p — > U pU^ . Henceforth, I 
will also assume h — 1. 

If the unitary operator U n is weak, that is, close to the identity, one can always find a 
Hamiltonian operator H n such that 

U n = exp {-iHJt} « 1 - iHJt . (10) 

Thus, one can easily recover the Schrodinger equation from a description in terms of unitary 
operators, 

5\^ n ) = |^n) - |^n-l) = (U n ~ l)|Vn-l> = -iif n |^ n _i)<St • (H) 



Exercise 1. Show that the most general unitary operator on a single q-bit can be written 

U — 1 cos ip + in • a sin ip — exp{iipn ■ a} , (12) 



modulo an irrelevant overall phase, where n = {n x ,n y ,n z } is a three-vector with unit norm 
and a = {a x , cr y , a z }. In the Bloch sphere picture this is a rotation by an angle if about the 
axis defined by n, with three real parameters (four if one includes the global phase). 

For two q-bits there is unfortunately no such simple visualization as the Bloch sphere. 
However, to specify any unitary transformation it suffices to give its effect on a complete set 
of basis vectors. I will consider only a fairly limited set of two-bit transformations in this 
paper, and no transformations involving more than two q-bits, but the simple formalism I 
derive readily generalizes to higher-dimensional systems. 

Let us examine a couple of examples of two-bit transformations. The controlled-NOT 
gate (or CNOT) is widely used in quantum computation; applied to the tensor-product basis 
vectors it gives 



tWr|00> 
£^cnot|01) 

t>CNOT|10> 

^cnot|11) 



|00> 
|01) 

111) 

|10) 



(13) 



or in matrix form 



CNOT 



/I 








°^ 





1 

















1 


Vo 





1 


o) 



(14) 



If the first bit is in state |0) this gate leaves the second bit unchanged; if the first bit is 
in state |1) the second bit is flipped |0) <-> |1). Hence the name: whether a NOT gate is 
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performed on the second bit is controlled by the first bit. In terms of single-bit operators 
this is f> CNO T = |0)(0| <g> i + <g> a x . 

Another important gate in quantum computation is the SWAP; applied to the tensor- 
product basis vectors it gives 

f^SWAP 
^SWAP 
^SWAP 
f^SWAP 

which in matrix form is just 



t^SWAP 



r\r\\ 

00) 


|aa\ 

= |oo) 


5 


01) 


= |10> 


) 


10) 


= |01> 




111) 


= 111) 


1 


/I 





o\ 





1 








1 


1 


Vo 





1/ 



(15) 



(16) 



As the name suggests, the SWAP gate just exchanges the states of the two bits: C/swapOV') <8> 

10)) = |0)®|^). 

CNOT and SWAP are examples of two-bit quantum gates. Such gates are of tremendous 
importance in the theory of quantum computation. These two gates in particular have an 
additional useful property. Note that the operator Uqnot — ^cnot> i- e -> ^ * s both unitary 
and Hermitian. This is also true of ?7swap- This is only possible if all of the operator's 
eigenvalues are 1 or —1. This also means that these operators are their own inverses, 
^cnot = ^swap = 1- (Among single-bit operators, the Pauli matrices & x ,y,z also have this 
property, as does any operator n • a.) 

This property of being both unitary and Hermitian makes it possible to define one- 
parameter families of two-bit unitary transformations: 



and similarly 



U C not(9) = exp [-iOUcNor} 

= 1 COS 6 — iUcNOT sm ® ) 

^swap(#) = exp [-iOUswAp] 

= i cos 9 — iUswAP sin 9 . 



(17) 



(18) 



Exercise 2. Show that any operator U which is both unitary and Hermitian satisfies 
exp(— i9U) = i cos 9 — ill sin 9. 

These families range from the identity 1 for 9 = to the full CNOT (or SWAP) gate for 
9 = 7r/2, up to a global phase. For small values S < 1, this is a weak interaction, which 
leaves the state only slightly altered. 

These families do have an undesirable characteristic, however. Suppose we apply 
^cnot(^) to a two-bit state of the form (a\0) + (3\l)) ® The new state is 
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#CNo T (0)(a|O> + 0\1)) ® |V) = (cwr*|0) + /3cos0|l» ® |V> 

sin6»|l> <g> (ff x |^)) . (19) 

This relative phase e je between |0) and |1) is a complication in the calculations that will 
follow. To avoid this problem, instead of ?7cnot(#) and £7swap(#) we will use interactions 
of the form Z a (0)U C n O t(0), where Z A (9) = (exp(—ia z 9/2) <g> 1), and similarly for £7swap- 
The single-bit rotation exactly undoes the extra relative phase produced by £/cnot(#), while 
changing nothing else. 

Before leaving the topic of two-bit gates, I should point out another representation which 
is sometimes useful. Any operator on two q-bits can formally be written 



6 = R iA 
i,j=0 



'31 



(20) 



where cq = 1, <3"i,2,3 = &x,y,z- The matrix elements can be calculated using the identity 



R ij = Tr{6(a i ®a j )} . 



(21) 



If O — is Hermitian, the matrix elements Rij will be real. 
If we write the CNOT in this form, it gives matrix elements 



R 



CNOT 



/ 1/2 


1/2 































V-l/2 


1/2 





0/ 



(22) 



Similarly, if we write the SWAP gate this way it has matrix elements 



R 



SWAP — 



/1/2 








^ 





1/2 














1/2 





V o 








1/2 J 



(23) 



III. STRONG AND WEAK MEASUREMENTS 

A. Projective measurements 

In the standard description of quantum mechanics, observables are identified with Hermi- 
tian operators O = OK A measurement returns an eigenvalue o n of O and leaves the system 
in the corresponding eigenstate \<f> n ), 0\(f> n ) = o n \(f> n ), with probability p n = \(<j) n \ilj)\ 2 . If a 
particular eigenvalue is degenerate, one instead uses the projector V n onto the eigenspace 
with eigenvalue o n ; the probability of the measurement outcome is then p n = (ip\V n \ip) and 
the system is left in the state V n \ip) / ' ^fv~n- F° r a mixed state p the probability of outcome n 
is p n = Tr{P n p} and the state after the measurement is V n pV n /p n . 
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Because two observables with the same eigenspaces are completely equivalent to each 
other (as far as measurement probabilities and outcomes are concerned), we will not worry 
about the exact choice of Hermitian operator O; instead, we will choose a complete set of 
orthogonal projections {P n } which represent the possible measurement outcomes. These 
satisfy 

^PnPn' = T^n^nn' 1 ^ ] 'Pn = 1 • (24) 
n 

A set of projection operators which obey (|24l) is often referred to as an orthogonal decompo- 
sition of the identity. For a single q-bit, the only nontrivial measurements have exactly two 
outcomes, which we label + and — , with probabilities p + and p_ and associated projectors 
of the form 

P± = = llkW±l 5 (25) 

this is equivalent to choosing an axis n on the Bloch sphere and projecting the state onto 
one of the two opposite points. All such operators are projections onto pure states. The two 
projectors sum to the identity operator, = 1. The average information obtained from 

a projective measurement on a q-bit is just the Shannon entropy for the two measurement 
outcomes, 

SWas = ~P+ l0g 2 P+ - P- l0g 2 p_ . (26) 

The maximum information gain is precisely one bit, when p + = p_ = 1/2, and the minimum 
is zero bits when either p + or p_ is 0. After the measurement, the state is left in an eigenstate 
of V±, so repeating the measurement will result in the same outcome. This repeatability is 
one of the most important features of projective measurements. 



B. Positive operator valued and weak measurements 

While projective measurements are the most familiar from introductory quantum me- 
chanics, there is a more general notion of measurement, the positive operator valued mea- 
surement (POVM)El. Instead of giving a set of projectors which sum to the identity, we give 
a set of positive operators E n which sum to the identity: 

(27) 

n 

The probability of outcome n is p n = (ip\E n \ip), or for a mixed state p n = Tr{E n p}. Unlike 
a projective measurement, knowing the operators E n is not sufficient to determine the state 
of the system after measurement. One must further know a set of operators A nk such that 

E n = Y,A ] nk A nk . (28) 

k 

After measurement outcome n the state is 
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p _> p > = (i/ Pn ) A nkP Al k 

k 



(29) 



This measurement will not preserve the purity of states, in general, unless there is only a 
single A n k for each E n (that is, E n = A^A n ), in which case = A n \ip) / -Jp^,. 

Because the positive operators E n need not be projectors, one is not limited to only 
two possible outcomes; indeed, there can be an unlimited number of possible outcomes. 
However, if a POVM is repeated, the same result will not necessarily be obtained the second 
time. Projective measurements are clearly a special case of POVMs in which results are 
repeatable. Most actual experiments do not correspond to projective measurements, but are 
described by some more general POVM. 

Conversely, it is easy to show that any POVM can be performed by allowing the system 
which is to be measured to interact with an additional system, or ancilla, and then doing 
a projective measurement on the ancilla. In this viewpoint, only some of the information 
obtained comes from the system; part can also come from the ancilla, which injects extra 
randomness into the result. This mathematical identification should not be pushed too far, 
however; in practice, trying to interpret a POVM as a projective measurement on an ancilla 
often adds complexity without improving understanding. 

A particularly interesting kind of POVM for the present purpose is the weak 
measurement. This is a measurement which gives very little information about the system 
on average, but also disturbs the state very little. 

Loosely speaking, there are two ways a measurement can be considered weak. Suppose 
we have a q-bit in a state of form ([!]), and we perform a POVM with the following two 
operators: 

£ = |0><0| + (l-6)|l><li=i 2 , 
E 1 = e\l)(l\=Al, 

a = |o)(o| + v / r^|i)(i| , 

A 1 = V~e\l)(l\, (30) 

where e <C 1. Clearly Eq and E\ are positive and Eq + E\ = 1, so this is a POVM. The 
probability po = (ip\E \ip) = 1— e\/3\ 2 of outcome is close to 1, while p\ = (i/j\Ei\ip) = e|/3| 2 is 
very unlikely. Thus, most such measurements will give outcome 0, and very little information 
is obtained about the system. 

Exercise 3. Show that the average information gain from this measurement (in bits) is 

S meas « e\f3\ 2 (l/ In 2 - log 2 e\(3\ 2 ) < 1 . (31) 

The state changes very slightly after a measurement outcome 0, with |0) becoming slightly 
more likely relative to |1). The new state is 

|-0o> = Ao^)/Vp~o = HO) + pVT^\l))/y/to . (32) 

However, after an outcome of 1, the state can change dramatically: a measurement 
outcome of 1 leaves the q-bit in the state |1). We see that this type of measurement is 
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weak in that it usually gives little information, and has little effect on the state; but on rare 
occasions it can give a great deal of information and have a large effect. 
By contrast, consider the following positive operators: 



1 + e 



|0)(0|+(^)|l)(l|=^o 
1 + ^Ddhi? 




(33) 

These operators also constitute a POVM. Both Eq and E\ are close to 1/2, and so are almost 
equally likely for all states, po ~ p\ ~ 1/2; the information acquired is approximately one 
bit. But the state of the system is almost unchanged for both outcomes, with the new state 
being [ipj) = A^/y/pJ, 

\ipo) = i 1 = (aVT+7|0> + pVT^~e\l 

^ a(l + e\(3\ 2 )\0) + (3(1 - e\a\ 2 )\l) , 
= i 1 — (aVT^~e\0) + (3VTT~e\l 



1 - e( a 



^ a(l - e\(3\ 2 )\0) + (3(1 + e\a\ 2 )\l) } (34) 

for outcomes and 1, respectively. For this type of weak measurement, the measurement 
outcome is almost random, but does include a tiny amount of information about the state. 
By performing repeated weak measurements and examining the statistics of the results, one 
can in effect perform a strong measurement; for the particular case considered here, the 
state will tend to drift towards either |0) or |1) with overall probabilities \a\ 2 and \(3\ 2 . 

We will see below specific instances of how a POVM can arise by the interaction of the 
measured system with an ancilla, and the subsequent projective measurement of the ancilla. 
Quantum trajectories are closely related to such indirect measurement schemes. 



IV. SYSTEM AND ENVIRONMENT 

In discussing quantum evolution it is usually assumed that the quantum system is very 
well isolated from the rest of the world. This is a useful idealization, but it is rarely realized 
in practice, even in the laboratory. In fact, most systems interact at least weakly with 
external degrees of freedom0lll. 

One way of taking this into account is to include a model of these external degrees of 
freedom in our description. Let us assume that in addition to the system in state in 
Hilbert space Tis with Hamiltonian Hs, there is an external environment in state \E) e "He 
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with Hamiltonian H E . The joint Hilbert space of the system and environment is the tensor 
product space Tis <%>Ti E , and the full Hamiltonian can be written 

H = H s ®i E + l s ®H E + H J , (35) 

where Hj is a Hamiltonian giving the interaction between the two. The system and environ- 
ment evolve according to the Schrodinger equation (|?p with this full Hamiltonian. Because 
of the presence of the Hj term, in general even initial product states \ip) ® \E) will evolve 
into states which cannot be written as product states: 

M®|tf)-H*)=£lV'i>®|tfi)- (36) 

3 

Such states are called entangled. If the system and environment are entangled, it is no longer 
possible to attribute a unique pure state to the system (or environment) alone. The best 
that can be done is to describe the system as being in a mixed state p sys , which is the partial 
trace of the joint state |^) over the environment degrees of freedom: 

Ps ys = Tr env {|^)(^|} 

= \^3)^j'\{ e k\Ej){E r \e k ) 

= Y l \^)^A(Ef\E j ) , (37) 

3,j' 

where the {|efc)} are some complete set of orthonormal basis vectors on the environment 
Hilbert space TC E . 

The entanglement of a state can be quantified. The most widely used measure for pure 
states is the entropy of entanglement. For a system and environment in a joint pure state, 
this is 

S E = -Tr{p s log 2 p 5 } , (38) 

which is the von Neumann entropy of the reduced density operator p$. This entropy is 
maximized by the density operator ps = ls/Tr sys {l,s}; such a density operator is called 
maximally mixed, and the joint pure state which gives rise to it is called maximally entangled. 

Exercise 4. Show that the von Neumann entropy is maximized for p = 1/Tr{i}. 

Systems which interact with their environments are said to be open. Most real physical 
environments are extremely complicated, and the interactions between systems and envi- 
ronments are often poorly understood. In analyzing open systems, one often makes the 
approximation of assuming a simple, analytically solvable form for the environment degrees 
of freedom. 

For the models of this paper, I will assume that both the system and the environment 
consist solely of q-bits. I will also assume a simple form of interaction, namely that the 
system q-bit interacts with one environment q-bit at a time, and that after interacting they 
never come into contact again; and that the environment q-bits have no Hamiltonian of their 
own, H E = 0. This environment may seem ridiculously oversimplified, but in fact it suffices 
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to demonstrate most of the physics exhibited by much more realistic descriptions. And in 
fact, the description in terms of q-bits is not a bad schematic picture of many environments 
at low energies. I summarize one such possible environment in Figure 2. 

For this type of model, the Hilbert space of the system is Tis — ^2 and the Hilbert space 
of the environment is He = <8> ^2 (£>•••. I will in general assume that all the environment 
q-bits start in some pure initial state, usually |0); however, in later elaborations of the model 
I will relax this assumption to include both other pure states and mixed-state environments, 
which resemble heat baths of finite temperature. 



V. THE MASTER EQUATION 

A system interacting with an environment usually cannot be described in terms of a pure 
state. Similarly, it is not usually possible to give a simple time-evolution for the system alone, 
without reference to the state of the environment. However, under certain circumstances 
one can find an effective evolution equation for the system alone, if the interaction between 
the system and environment has a simple form and the initial state of the environment has 
certain properties. 

In cases where such an effective evolution equation exists, it will not be of the form ([7j). 
Instead, the evolution equation will be completely positive, which is a weaker condition than 
unitarity. It must take density matrices to density matrices, but not necessarily pure states 
to pure states. To be tractable, such an equation must also generally be Markovian: it must 
give the evolution of the density operator solely in terms of its state at the present time, 
with no retarded terms. (While it is possible to find retarded descriptions for more general 
non-Markovian open systems, thev tend to be difficult to solve, even numerically.) The most 
general such Markovian equation^ is a master equation in Lindblad form 

^ = -i[^,p]+E[WI-(V2)^*p-(l/2)p^*] , (39) 

k 

where H$ is the Hamiltonian for the system alone and the Lindblad operators give the 
effects of interaction with the environment. 

What about the kind of discrete time evolution we've been assuming? Suppose that a 
system and environment in an initial product state undergo some unitary transformation U : 

|tf ) = |V) <8> \E) = U (\if>) ® \E)) . (40) 

How does the state of the system alone change? We can write any operator on Tis <8> 'He as 
a sum of product operators 

U^AjQBj. (41) 
3 

The state of the system after the unitary transformation is 

f/ s = Tr m {U\*)(*\fr} 

= £Tr env {(i,#)(Vl4) ® (Bj\E)(E\By)} 
3,3' 

= Y,A^)^\A] I (E\B},B J \E) . (42) 

3,3' 
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The self-adjoint matrix Mjj> = (E\Bj,Bj\E) has a set of orthonormal eigenvectors jl k = {p k j} 
with real eigenvalues X k such that 

E Mjj>n kj > = \ k /j,kj Mjf = E ^kVkjUlf ■ (43) 
j> k 

We define new operators O k by 



O k = \J\kY, PkjAj . (44) 
j 

Exercise 5. Show that the unitarity of U implies 

E 0\O k = Is . (45) 
k 

The above expression then simplifies to 

p' s = Tr env {?7|*><*|C/t} = J2d k \^)^\dl . (46) 

k 

This expression is quite like the outcome of a POVM, but without determining a particular 
measurement result; the system is left in a mixture of all possible outcomes. 

Suppose now that the system interacts with a second environment q-bit in the same way, 
with the second bit in the same initial state \E), and that the Hilbert space of the of two 
environment bits is traced out. The system state will become 

p" s = £ d k2 d kl mWdlo~l ■ (4T) 

fcl,£l2 

After interacting successively with n environment bits the system state is 

pP= E O kn ---d k M)^\Ol---Ot n . (48) 

k i , . . . , hn 

This evolution is a type of discrete master equation. We see clearly that an initially pure 
state will in general evolve into a mixed state. Depending on the O k , it may or may not 
converge to a unique final state. 

Let's consider the special case when the unitary operator U is close to the identity, 

U = exp j-ieE^j ® Bj\ , (49) 

where e C 1, Tr-fA-A,-} ~ Tr{BjBj} ~ 0(1), and the A,- and Bj are Hermitian. Expanding 
to second order in e, we see that the new density matrix for the system is 
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P ' s = Tr env {U\*)(*\tf} 

*\^)^\-zeY,lA v \^(ME\B 3 \E) 

j 

+^E( s l4'4l^> 

33' 

x (2A^){^\A r - A f A^)(^\ - \1>)U>\Aj,Aj) . (50) 
Let's make the simplifying assumption that the first-order term vanishes: 

J2A ] (E\B,\E) = 0. (51) 

3 

Later in the paper we supplement ?7cnot(#) and £/swap(#) with one-bit gates on the system 
to remove this first-order term, which simplifies the evolution. 

Let the duration of the interaction be St. We can define a matrix Mjf 

M jf = {E\B r B 3 \E) , (52) 

with orthonormal eigenvectors jlj- = {/%} and real eigenvalues just as in ([43]), and define 
operators 



St 



Y^VkjAj, (53) 



similarly to (44]) above. In terms of these operators, equation ( ]50|) can be re-written 

^ = £ [LkpLl - {l/2)LlL k p - {l/2)pL\L k ] , (54) 

which has exactly the same form as ([39]) with a vanishing Hamiltonian. If the first order 
term of flBTp does not vanish, equation (Bl) can pick up an effective system Hamiltonian, 
though some care must be exercised about the size of the various terms and their order 
in St. But for the simple case we are considering, we can recover a continuous evolution 
equation as a limit of successive brief, infinitesimal transformations, just as we recovered the 
Schrodinger equation from a succession of weak unitary transformations. 

If we wanted to express the operators Ok of equation (H) in terms of the L k for this 
infinitesimal transformation, then to second order in e we would find Ok = Lk\fSi, plus 
one additional operator O = 1 — (1/2) J2k L\LkSt. Substituting these operators in 
reproduces ([54]). 



VI. ENVIRONMENTAL MEASUREMENTS AND 
CONDITIONAL EVOLUTION 

Suppose that the system q-bit is in initial state a|0) + /3|1) and it interacts with an 
environment q-bit in state |0), so that a CNOT is performed from the system onto the 
environment bit. The new joint state of the two q-bits is 
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\V') = a\0) s ®\0) E + (3\l)s®\l)E- (55) 
If we trace out the environment bit the system is in the mixed state 

p 5 = |a| 2 |0)(0| + |/?| 2 |l)(l|. (56) 

Now suppose that we make a projective measurement on the environment bit in the usual 
z basis. With probability po = \a\ 2 (pi = |/3| 2 ) we find the bit in state |0) (|1)), in which 
case the system bit will also be projected into state |0) (|1)). If the system interacts in the 
same way with more environment bits, which are also subsequently measured in the z basis, 
the same result will be obtained each time with probability 1. This scheme acts just like a 
projective measurement on the system. 

This need not be the case in general. Suppose that instead of a CNOT, a SWAP gate is 
performed on the two q-bits. In this case the new joint state after the interaction is 

m = \0) s ®(a\0) E + (3\l) E ) ■ (57) 
A subsequent z measurement on the environment yields or 1 with the same probabilities 
as before, but now in both cases the system is left in state |0). Further interactions and 
measurements will always produce the result 0. Clearly, the choice of measurement on 
the environment bit makes no difference to the state of the system — it will be |0) for any 
measurement result, or for no measurement at all. This is rather like spontaneous decay 
of a system which is initially in a superposition of the ground and excited states. If we 
think of |1) as an excited state, the system passes on an excitation to the environment with 
probability \(3\ 2 and is afterwards left in the ground state. Note that Ucnot can produce 
entanglement between two initially unentangled q-bits, while {/swap cannot. 

In the case of the CNOT interaction, what happens if we vary our choice of measurement? 
Suppose that instead of measuring the environment bit in the z basis, we measure it in the 
x basis \x±) = (|0) ± |l))/\/2? In terms of this basis we can rewrite the joint state after the 
interaction as 

I*') = ^(«|0> + ® \x+) + ^(a|0> - (3\1)) ® l*-> • (58) 

The + and — results are equally likely; the + result leaves the system state unchanged, while 
the — flips the relative sign between the |0) and |1) terms. Aside from knowing whether a flip 
has occurred, this measurement result yields no information about the system state; if this 
state is initially unknown, it will remain unknown after the interaction and measurement. 
By changing the choice of measurement one can go from learning the exact state of the 
system to learning nothing about it whatsoever. 

One important thing to note is that, whatever measurement is performed and whatever 
outcome is obtained, the system and environment bit are afterwards left in a product state. 
No entanglement remains. Thus, the environment bit can be further manipulated in any 
way, or discarded completely, without any further effect on the state of the system. Because 
of this, it is unnecessary to keep track of the environment's state to describe the system. 
This is a somewhat unrealistic approximation, arising from the assumption of projective 
measurements — in real experiments, some residual entanglement with the environment al- 
most certainly exists which the experimenter is unable to measure or control. It does, 
however, enormously simplify the description of the system's evolution. I will later relax 
this assumption, but for now I make use of it. 
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VII. STOCHASTIC SCHRODINGER EQUATIONS WITH 

JUMPS 



Far more complicated dynamics result if we replace the strong two-bit gates of the pre- 
vious section with weak interactions, such as {7cnot(0) or i7swAp(0) for 6 <C 1. In this case 
the measurements on the environment bits correspond to weak measurements on the system, 
and the system state evolves unpredictably according to the outcome of the measurements. 
Weak interaction with the environment is the norm for most laboratory systems currently 
studied; indeed, it is only the weakness of this interaction which justifies the division into 
system and environment in the first place, and enables one to treat the system as isolated 
to first approximation, with the effects of the environment as perturbations. 

This problem is simple enough that we can analyze it for a general weak interaction. 
Suppose that the system and environment q-bits are initially in the state |\&) = \ip) s ® |0) E . 
Let the interaction be a two q-bit unitary transformation of the form 

C/ = exp r <g> £,j , (59) 

where ^ < 1, and Tr{A,} ~ Tr{Bj} ~ 0(1), so we can expand the exponential in powers 
of 9. The new state of the system and environment is 

U\V) = |V) s ® \0)e ~ id E A M)s ® B 3 \0) E 

j 

- e lY, A r^W)s®B j ,B j \Q) E + 0{e") . (60) 

33' 

Just as in section 5, we define the matrix Mjj> = (0\Bj>Bj\0) with eigenvectors /4 and 
real eigenvalues and use it to find a new set of operators 



L k = \j 6 -^Y. VkjAj , (61) 

where St is the time between interactions with successive environment q-bits. Since we can 
decompose the matrix Mjj/ in terms of two vectors 

M 3f = (0|4-,]0)<0|4|0) + (0|4-,|1><1|4-|0) = u j/Uj + v} Vj , (62) 

it has at most two nonvanishing eigenvalues. We can make the additional simplifying as- 
sumption, as in section 5, that 

£i,(0|4|0> = 0. (63) 

3 

This leaves only a single nonvanishing Lindblad operator, 

£ = y^£i,(ii4io>, (64) 



St 
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and no effective Hamiltonian term. In terms of L, the joint state of the system and environ- 
ment is 

U\V) = (i - (l/2)ti L6t)\i/>) s ® |0) B - iV6iL\^j) s ® \1) E + 0{6 3 ) . (65) 



The time dependence of L ~ 6/ySt may seem a bit strange at first. It is chosen so 
that the master equation ([54]) properly approximates a time derivative. This does, however, 
imply a somewhat odd scaling. The strength of the interaction between the system and 
environment q-bits can be parametrized by 6 /St. This implies that if the two bits interacted 
for only half as long, the Lindblad operators would have to be redefined, L — > L/y/2. In the 
limit as St — > 0, the right-hand-side of fl5"4|) would vanish, and the system would not evolve. 
This is actually true, and is an example of the so-called Quantum Zeno Effect^. Most such 
master equations have some microscopic timescale built in — for instance, a characteristic 
collision time — and are not really valid at shorter timescales. These equations are used to 
describe the system at times long compared to this short scale, which in the case of our 
model would be after the system had interacted with a number of environment q-bits. 

Once the interaction has taken place, we measure the environment q-bit in the {|0), |1)} 
basis. Using the assumption fl63|), the probabilities p ,i of the outcomes are 



Po = 1 - 2 £ (M-'^#><oi4-'|i><il4-|o> 
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= 1 - (ip\P L\ip)St 

= 1-Pi- (66) 
The change in the system state after a measurement outcome on the environment is 

sty) = |v/> - m = -i {pi - ipLj) Wt . (67) 

Since the result is highly probable, as the system interacts with successive environment 
bits in initial state |0) £ , most of the time the system state will evolve according to the 
nonlinear (and nonunitary) continuous equation (|6"7D. 

Occasionally, however, a measurement result of 1 will be obtained. In this case, the state 
of the system can change dramatically: 

W - W) = ^H= • (68) 



Because these changes are large but rare, they are usually referred to as quantum jumps. 

These two different evolutions — continuous and deterministic (after a measurement result 
0) or discontinuous and random (after a measurement result 1) — can be combined into a 
single stochastic Schrddinger equation: 

6\4>) = M-\4>) = -l(tiL-(ftL))\il>)6t 

-i)\i>)SN, (69) 
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where SN is a stochastic variable which is usually 0, but has a probability pi = (if)\UL\ifj)St 
of being 1 during a given timestep St. The values of SN obviously represent measurement 
outcomes; when the environment is measured to be in state |0) the variable is SN = 0; 
when the measurement outcome is 1, SN = 1. The equation flSPp combines the two kinds 
of system evolution into a single equation. Most of the time the SN term vanishes, and 
the system state evolves according to the deterministic nonlinear equation flo"7|); however, 
when SN = 1 the SN terms completely dominate, and the system state changes abruptly 
according to (|68l). Generically, (|69|) can also include a term for Hamiltonian evolution, but 
this is eliminated by the assumption (|63|). 

We can summarize the behavior of the stochastic variable by giving equations for its 
ensemble mean, 

(SN) 2 = SN , M[SN] = (L^L)St , (70) 

where M[] represents the ensemble mean over all possible measurement outcomes, and () is 
the quantum expectation value in the state l^). 

This analysis is rather abstract. Let us see how this works with a particular choice of 
two-bit interaction. Suppose that the system interacts with a succession of environment q- 
bits, via the transformation Z s (9)Ucnot{9), where Z s {9) = exp(— i9a z /2)s <8> 1e an d # <C 1; 
and that after each interaction the environment bits are measured in the z basis. If the 
environment bits all begin in state |0), and the system in state ([]]), after interacting with 
the first environment q-bit the two bits will be in the joint state 

\W) = a |oo) + /3cos0|lO) - i/3 sin 0| 11) , (71) 

modulo an overall phase. This state is entangled, with entropy of entanglement Se ~ 
— \a(3\ 2 9 2 log 2 \(af3\ 2 9 2 ). If the environment is then measured in the z basis, it will be found 
in state |0) with probability po = \a\ 2 + |/3| 2 cos 2 6* ~ 1 — \(3\ 2 9 2 and in state |1) only with 
the very small probability p\ = \(3\ 2 6 2 . After a result of 0, the system will then be in the 
new state 

|V>o) = (a|O)+/3cos0|l))/v^ 

w a(l + \P\ 2 9 2 /2)\0) + p(l - \a\ 2 9 2 /2)\1) , (72) 

where we have kept terms up to second order in 9. The amplitude of |0) increases relative 
to |1). On the other hand, if the outcome 1 occurs, the system will be left in state |1). 
Further interactions and measurements will not alter this. We see that this combination of 
interactions and measurements is exactly equivalent to a weak measurement on the system, 
of the first type discussed in section 3.2. This bears out the fact that all POVMs are 
equivalent to an interaction with an ancillary system followed by a projective measurement, 
the well-known "Neumark's Theorem." 

Suppose that the system interacts successively with n environment q-bits, each of which 
is afterwards measured to be in state |0). The state of the system is \if) n ) = a n \0) + /3 n \l), 
where 

-=(l-2-J -«e*p(-n*/2)-, (73) 
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which, together with the normalization condition \a n \ 2 + \f3 n \ 2 = 1, implies 

12 



\(X 



|2 



\Ct\ 



\a\ 2 + \(3\ 2 exp(-n6 2 ) ' 
2 |/?| 2 exp(-^ 2 ) 

|a| 2 + \{3\ 2 exp(-n6 2 ) ' 1 J 

We see that \a n \ 2 — > 1, which is what one would expect, since after many measurement 
outcomes without a single 1 result, one would estimate that the |1) component of the state 
must be very small. The probability of measuring at the nth step, conditioned on having 
seen Os at all previous steps, is po(n) ~ 1 — \(3 n \ 2 9 2 — > 1. 

Exercise 6. Show that the the probability of getting the result at every step as n — > oo, 
which leads to the system asymptotically evolving into the state |0), is 

oo 

f[po(n)n\a\ 2 . (75) 

n=0 



The probability of at some point getting a result 1 and leaving the system in state |1) is 
therefore 1 — |a| 2 =|/3| 2 . The effect, after many weak measurements, is exactly the same as 
a single strong measurement in the z basis, with exactly the same probabilities. 

This evolution is described exactly by the equation ( |69|) with Lindblad operator 

L=4¥jbt\\){\\. (76) 
In this case, the equation simplifies to 

W) - W) = -y " I W)| 2 ) m + - i) , (77) 

where M[5N] = 6 2 \{l\ijj)\ 2 = 6 2 \j3\ 2 = p±. Equation ( |69D may seem unnecessarily compli- 
cated. After all, L is Hermitian, so the distinction between L and D is unimportant; also, 
L is proportional to the projector |1)(1|, which simplifies (|69|) considerably. However, many 
other systems exist in which L is not so simple, and which still obey equation (|69|). 

We can readily find such an alternative system. Suppose that t/swAp(69 is used instead of 
^cnot(^) in the two-bit interaction. The system and environment bits after the interaction 
are in the state 

= a |00) + /ScosflllO) - z/3sinfl|01) . (78) 

(Note that unlike the strong swap gate ?7swap, the weak £/swap(#) can produce entangle- 
ment.) From (|78|) we see that just as in the CNOT measurement of on the 
environment leaves the system in a slightly altered state — in fact, exactly the same state 
that is produced by a result in the CNOT case. The probabilities po and p± for the two 
outcomes are exactly identical to those for the CNOT case, po = 1 — p\ ~ 1 — |/3| 2 # 2 . How- 
ever, unlike the CNOT case, a measurement outcome 1 will put the system into the state 
|0) rather than |1). 
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This process is quite like spontaneous decay. Suppose that |0) is the ground state and |1) 
is the excited state. Initially, the system is in a superposition of these two energy levels. Each 
timestep, there is a chance for the system to emit a quantum of energy into the environment. 
If it does, the system drops immediately into the ground state |0), while the environment 
goes into the excited state |1); if not, we revise our estimate of the system state, making it 
more probable that the system is unexcited. This is exactly similar to quantum jumps in 
quantum optics, in which a photodetector outside a high-Q cavity clicks if a photon escapes, 
and gives no output if it sees no photoni'lfi. A measurement result of 1 on the environment 
corresponds to the photodetector click. 

This model obeys exactly the same stochastic Schrodinger equation fl69|), but with the 



Lindblad operator L = ^9 2 /St\0)(l\. Because this operator is not Hermitian, the distinction 

between U and L is important in this case. We compare these two stochastic evolutions in 
Fig. 3a and 3b. In both cases, the coefficient \(3\ 2 decays steadily away; however, in some 
trajectories the state suddenly jumps either to |1) or |0), while others continue to decay 
smoothly towards |0). 

What happens if we don't actually measure the environment bits? In this case the 
system q-bit evolves into a mixed state which is the average over all possible measurement 
outcomes. If we denote by £ the n outcomes after n steps of the evolution, the probability 
of those outcomes, and |^) the state conditioned on those outcomes (i.e., the solution of 
the appropriate stochastic Schrodinger equation), then the system state will be a density 
operator 

P (n) = £I^M£X^I- (79) 

If we start with a state p, allow it to interact with exactly one environment q-bit, and average 
over the possible outcomes and 1 with their correct probabilities, we will find that the new 
state p' is 

p' = j2 A kP A\ , (80) 

k=0 

where for ?7cnot 

A = |0)(0| + cos0|l)(l| , 

Ai = sin0|l)(l| . (81) 
In the limit of small 6^1 this evolution gives an approximate Lindblad master equation 

^ = L P D - l -VL P - \ p Dl . (82) 

The Lindblad operator L is the same as in the jump equation for the CNOT interaction, 
([761) ■ Similarly, for the SWAP interaction, the evolution equations are identical, but with 



A 1 = sin0|O)(l| and L = J6 2 /5t\0)(l 
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We should note that unlike the stochastic trajectory evolution, the evolution of the 
density matrix is perfectly deterministic, with no sign of jumps or other discontinuities. 
This is generally true; the jumps appear only if measurements are performed. 

Believe it or not, the models analyzed in this section have already given almost all 
essential properties of quantum trajectories. Quantum trajectory equations arise when a 
system interacts with an environment which is subsequently measured; the equations give 
the evolution of the system state conditioned on the measurement outcome. They take 
the form of stochastic nonlinear differential equations; the stochasticity arises due to the 
randomness of the measurement outcomes, and the nonlinearity due to renormalization of 
the state. Averaging over all possible measurement outcomes recovers the deterministic 
evolution of the system density operator. The next two sections examine elaborations of 
this basic scheme, and section 10 discusses how these results fit into the foundations of 
quantum mechanics. 

VIII. STOCHASTIC SCHRODINGER EQUATIONS WITH 

DIFFUSION 

In discussing weak measurements we noted that some usually produced only minor ef- 
fects, but occasionally caused a large change in the state, while others changed the state 
little regardless of the outcome. The former case is exactly that of trajectories with jumps: 
the state usually changes slowly and continuously, but occasionally makes a large jump. 
The latter case, then, should correspond to states which evolve by small but unpredictable 
changes — in other words, by diffusion rather than jumps. 



A. Random unitary diffusion 

We can easily demonstrate such evolution using the model we've already constructed. 
Suppose that the system begins in a state (|IJ), and interacts with a succession of environment 
bits initially in state |0) by the unitary transformation Z s (9)Ucnot{9) for 9 <C 1. However, 
rather than measuring the environment bits in the z basis, we measure them in the x basis, 

k + ) = (|0) + |l))/v / 2, |x_> = (|0>-|l»/v^. (83) 
In terms of this basis, after interaction the system and environment bit are in the state 

|*') = -L ( a |0> + /3e"^|l>) |x + > + -L ( a | ) + /fe+"|l>) . (84) 

We see that the two measurement outcomes + and — are equally likely, and produce only a 
small effect on the system state: namely, if an outcome + occurs the relative phase between 
|1) and |0) is rotated by —9, and if a — occurs it is rotated by +9. This measurement 
scheme, then, has the effect of performing a random unitary transformation on the system. 
The stochastic Schrodinger equation for this model is 

S\tj)) = - |V>) = iL\ip)5W , (85) 
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where L = ^9 2 /St\l)(l\ is the same as in the earlier case of CNOT jumps, and SW is a 

stochastic variable taking the values SW = ±ySt with equal probability. We express this 
stochastic behavior by giving the ensemble means 

M[8W] = , SW 2 = St . (86) 

This particular normalization of the stochastic variable is useful because it scales properly 
with time. If, instead of using interactions with single environment q-bits in intervals St, we 
let the system interact with n bits over an interval At = nSt, then the equation (1851) remains 



completely unchanged, except for St being replaced by At in ( pq ) and (|86|). 

What if we average over these measurement outcomes? In this case, the density matrix 
evolution is still given by (BIT), with operators 



A = i=(|0><0|+e-* e |l><l 



i 1 = i=(|0)(0|+e^|l)(l|) . (87) 

This looks quite different from the evolution given by (|81|), but it is actually identical, as 
it must be: the same system state must be obtained by tracing out the environment bit, 
regardless of the basis in which we choose to carry out the trace. At the level of density 
matrices, the nonunitary quantum jumps and the unitary diffusion are exactly the same. 

Exercise 7. Show explicitly that the new density matrix 

p' = J2A kP Al (88) 

k 

obtained using the operators defined in fl37|) is the same as that using the operators defined 
in (BID. 



This is an example of a single master equation having different unravelings, that is, 
different stochastic evolutions QB5p and ( |8"5| ) which give the same evolution on average. Any 



set of measurements done on the environment must give the same average evolution of the 
system state, and therefore gives an unraveling of the same master equation. Of course, some 
unravelings may give simpler versions of the stochastic system state evolution than others. 
For some systems, it may be much easier to solve the stochastic evolution than the master 
equation; one can then average over many different stochastic evolutions to approximate 
the master equation evolution. This is called the Monte Carlo Wavefunction technique!. 
Even in cases where one has no direct access to information about the environment, it 
may be worthwhile to invent a fictitious set of measurements which would give a simple 
stochastic evolution, purely as a numerical technique to solve the master equation. In 
this case, however, one should not ascribe any fundamental significance to the individual 
solutions of the stochastic trajectory equations. 
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B. Nonunitary diffusion 



We can get diffusive behavior from this system in a different way. Suppose once again 
that the interaction is Zs(9)Ucnot(9) and the environment is measured in the z basis, but 
this time let the environment q-bits be initially in the state 

|y-) = (|0)-^|l))/v / 2. (89) 
In this case, the joint state after the interaction is 

\^') = 4=(>|0> + /3(cos0 - sin0)|l))|O) 

— ^=(a|0) +/3(cos0 + sin0)|l))|l) . (90) 
v2 

For small 9 <C 1 the two measurement outcomes are almost equally likely, po,i = (1/2) = F0|/3| 2 , 
and the conditioned states of the system are, to second order in 9, 

|^o)^«(l + |/3| 2 + 3|/3| 4 2 /2)|O) 

+(3(1 - \a\ 2 9 - \(3\ 2 9 2 - 9 2 /2 + 3|/3| 4 2 /2)|1) , 

|^)^«(1-|/?| 2 + 3|/?| 4 2 /2)|O) 

+/3(1 + |a| 2 - \(3\ 2 9 2 - 9 2 /2 + 3|/3| 4 2 /2)|1) . (91) 

This is essentially the same as the second weak measurement scheme (|33|) in section 3.2. 
If the system is initially in either state |0) or |1) it is unchanged, but all other states will 
diffuse along the Bloch sphere until they eventually reach either |0) or |1). Once the state 
is within a small neighborhood of |0) or |1) it is unlikely to diffuse away again. Because of 
this, this scheme can also be thought of as an indirect measurement of |0) vs. but in 
this case, the measurement result can only be determined by gathering a very large number 
of outcomes. If, after the system has interacted with many environment bits, we find that 
the number of 1 outcomes slightly exceeds the number of outcomes, we would conclude 
that the system is in state |1); if not, that the system is in state |0). 

This scheme is somewhat analogous to homodyne measurement in quantum optics@@. 
In homodyne detection, the weak signal from the system is coherently mixed with a strong 
local oscillator in a beam splitter; the two ports are fed into separate photodetectors. The 
output from each detector is almost entirely due to the local oscillator field; because this 
field is so much stronger than the signal from the system, almost all of the photons which 
arrive at the photodetectors come from the local oscillator, and very little information about 
the system is obtained in any one detector "click." To gain information about the system, 
one looks at the small difference between the photocurrents from the two detectors over the 
course of many "clicks." Similarly in our model, over the course of many events one can 
deduce information about the system state from the small excess of 1 results over results. 
The analogy to homodyne measurement is not as exact as the analogy between the quantum 
jump equation and spontaneous emission, but it is still suggestive. 

For this nonunitary diffusion, the stochastic Schrodinger equation is different than in 
the earlier quantum jump equation (|69|). Expanding the solution to second order in 9, the 
equation is 
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W)-\^) = \{?>\(L)\ 2 -2(D)L-Dl) \i;)5t 

+{L - (L))\ij)6W , (92) 

where once again L = J0 2 /5t\l){l\, and just as in the random unitary case of section 8.1 

the stochastic variable 5W takes the values ^fVfrt with probabilities po,i- Because in this 
case the measurement outcomes are not equally probable, this stochastic variable SW does 
not have vanishing mean: 

M[5W] = 2{D)5t , 6W 2 = 5t . (93) 

We can simplify this equation by replacing this stochastic variable by one that does have a 
vanishing mean (to second order in 9). Define 

5Z = 5W - M[8W] , 
M[5Z] = , 5Z 2 = St + 0(5t 3/2 ) . (94) 

Using this new stochastic variable the equation becomes 

m = w) - m « ({ti)L - l|(L)i 2 - \m 

+(L-(L)M5Z. (95) 

This is the quantum state diffusion equation with real noisM. Solutions of this equation are 
plotted in Figure 4. 

If we average over the measurement outcomes to give a mixed-state evolution for the 
system, the new density matrix is given by ( |30D with operators 



i o =^(|O><O| + (cos0-sin0)|l><l|) , 

h = -= (|0)(0| + (costf + sin0)|l)(l|) . (96) 
V2 



Exercise 8. Using the operators (pq), show that the density operator evolution is exactly 
the same as that given by operators ( |STD and (|STD, despite the fact that the environment's 
initial state is quite different in this case, and that the quantum trajectory evolution is also 
very different. We can see from this that the system evolution by itself is not sufficient to 
uniquely determine the nature of the environment. A particular master equation can result 
from many different environments and interactions. 



IX. INCOMPLETE INFORMATION AND STOCHASTIC 
MASTER EQUATIONS 

So far, we have only considered models in which both the system and environment are 
in a pure state, and measurements are done which leave that the case. The experimenter 
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has complete information about the state of the system and environment at all times. It 
is possible to generalize this in a numbers of ways. We may have only partial information 
about the state of the system, so that it is in an initial density matrix rather than a pure 
state. We may have only partial information about the environment; for instance, it might 
be in an initial thermal state. Or the measurements we perform on the environment may, 
themselves, give only partial information. In a realistic experiment, of course, all three of 
these may be true. We will consider each of them in turn. 



A. Trajectories for a mixed system state 

A pure state represents the maximum information that one can possess about a system. 
If a system is in a pure state, that implies that there is some measurement which could 
be performed on the system which would have a definite outcome. If we have only partial 
information about a system, we describe it by a mixed state. If a system is in a mixed state, 
there is no measurement which has a definite outcome. Systems can be mixed because we 
lack information about how they were prepared, or because they are entangled with other 
systems, or both. The von Neumann entropy is the most useful measure of how mixed a 
state is; it vanishes for a pure state, and is maximized by the maximally mixed state. 

For simplicity, let's suppose that we have no information at all about the state of the 
system; it is in the maximally mixed state p — 1/2, and has von Neumann entropy S = 1. We 
let it interact with an environment q-bit in some initial state, and measure the environment 
bit in some basis afterwards. 

For these measurement schemes, the state of the system becomes 

p -> p t = AipAl/pi , (97) 

with probability 

Pl = Tr{pijl t } . (98) 

If the interaction is Zs(9)Ucnot(@), as in the first quantum jump equation of section 7, the 
environment bits are initially in state |0), and are afterwards measured in the z basis, the 
operators (|97D are given by (|8lD . Just as for the pure state, an outcome 1 leaves the system 
in the state |1), while an outcome leaves the state only slightly changed: 

p ^i((l + ^ 2 /2)|0)(0| + (l-^/2)|l)(l|) , Po ^l-e 2 /2, 

Pl = |l)(l|, Pl ^ 6 2 /2. (99) 

In either case, the state is no longer completely mixed. Some information has been 
acquired about the system. We can assess how much information this is by calculating the 
average von Neumann entropy S of the state after the measurement: 

AS= 1-3= l + $>Tr{ Pi log 2 p 4 } , (100) 
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which in the case of our jump trajectories is AS ~ # 2 /(21n2). How much information does 
the measurement generate? This is given by the Shannon entropy (|26"D of the measurement, 



S m eas~-(#72)l0g 2 #72. (101) 

(Note that ini, Soklakov and Schack call S the average conditional entropy and denote it 
H, and call S meas the average preparation information and denote it 7.) 

We see that S meas > AS, which of course makes sense; we can never obtain more in- 
formation on average about the system than the measurement produces. If we look at the 
ratio of the two, though, 



S„ 



AS 



l-lnfl 2 /2, (102) 



we see that as 8 2 — > the information produced by the measurement is much, much higher 
than our gain in knowledge about the system. Thus, most of the information we receive is 
just random noise, telling us not about the state of the system but about the measurement 
process itself. 

If, instead of the jump description (^l|) given by measuring the environment in the z 
basis, we use the unitary diffusive description given by ( |3TD which arises from measuring 
the environment in the x basis, the situation is very different. Because Aq and A\ are in 
this case proportional to unitary operators, the initial maximally-mixed state p = 1/2 is 
left completely unchanged by this trajectory. Absolutely no information is gained about 
the system, AS = 0. This lack of information is not special to the maximally-mixed state. 
Because the two measurement outcomes are equally likely, po — pi — 1/2, each measurement 
result represents S'meas = 1 bit of information, but this tells one nothing about the state of 
the system; it is all random noise. 

Consider now the second diffusive case described in section 8.2, in which the environment 
bits begin in the state \y_) and, after interacting with the system, are measured in the z 
basis. In this case, with evolution operators given by (j96D, an initially maximally-mixed 
state is left in one of the states 



with probabilities 



Po,i^((l±0)|O)<O| + (lT0)|l)<l|) , (103) 



p o ,i~~(l±0). (104) 



Exercise 9. Show that, once again, most of the information in this measurement outcome 
is random noise, with the ratio AS/S meas — > as 9 — > 0. 



B. Trajectories for a mixed environment 

Suppose that, rather than the system being initially mixed, the environment bits are in 
a mixed state p^? For simplicity let's choose the diagonal state 
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Pe = w \0)(0\+w 1 \1)(1\ , (105) 

with wq + W\ = 1. Let an initially pure system in state (|J) interact with the environment 
bit via Zs{0)Uc7ssot{@)- The joint state of the system and environment after the interaction 
will be 

Ptot = |a| 2 w |00)(00| +a(3* Wq cos 9 1 00) (10 1 — iotj3*wo sin #|00) (11| 
+ |«|V|01)(01| -za/3*wisin0|Ol)(lO| +a(3*w 1 cosd|01)<ll| 
+a*/3w cos 0| 10) (00| + ia*/3wi sin 9\ 10) (01 1 
+ |/3| 2 (w cos 2 ^ + wi sin 2 6>)|10><10| 
+z|/?| 2 (wi - w ) cos 9 sin 0| 10) (1 1| + ia*/3w sin6>|l 1) (00| 
+a* / 9wiCos0|ll)(Ol| -i|/3| 2 (wi - w ) cos#sin#|ll)(10| 
+ |/3| 2 (wicos 2 ^ + wosin 2 e)|ll)(ll| . (106) 

For small 9 -C 1, if we measure the environment in the z basis the system will be left in 
one of the states 

p ,i « (1 -^l^lVoM,!)^')^'! + (^ 2 |^lVo/^o,i)|l)(l| , (107) 
with probabilities 

p 0ll ««)o 1 i±eWK-Wo) I (108) 



where 



=a(l + ^ 2 | / 9| 2 /2)|0)+/5(l-^ 2 |a| 2 /2)|l) . (109) 



These states po,i are both mixed; the uncertainty in the state of the environment has led to 
uncertainty in the state of the system. The system evolution is perturbed by noise from the 
environment. 

This evolution cannot be described by a stochastic Schrodinger equation; rather, one 
must use a stochastic master equation, which gives the evolution of a density matrix condi- 
tioned on the measurement outcomes. The master equation takes the form 



Ps ~ Ps = ~ 
+ 



l {DL-(DL), Ps }(^(l-F) + ^F)st 
2 1 J \Po Pi J 



Lp s D-{DL)ps] (^(l-F) + ^-F)5t. (110) 
J \Po Pi J 



Here {A, B} = Ab + bA is the anticommutator, and (O) = Tr{Op}. The stochastic variable 
F takes the value with probability po and 1 with probability p\. We see then that this 
equation ( |110|) reduces to (|82|) in the mean, and hence gives a partial unraveling of that 
equation in terms of mixed states. When w\ — > this master equation exhibits jump-like 
behavior; when Wi ~ w it is diffusive. If w = w\ = 1/2 then ( |110| ) reduces to the mean 



equation (32), since the measurements yield no information whatsoever about the system. 



(Note that this needn't be true for all possible interactions. For instance, if we use £/swap(0) 



27 



instead of C/cnot($)) it is possible to obtain information about the system even using a 
maximally mixed environment, as we see below.) 

If the system interacts with n environment bits, each of which is subsequently measured 
in the z basis giving a measurement record £, the system will be left in the state 

where depends on the exact sequence of measurement results, and \ip n ) is given by (|74|) 
and obeys \ip n ) — > |0) as n — ► oo. The state will tend to evolve towards either |0)(0| or 
|1)(1| unless w — w± — 1/2. So we see that even when a state becomes mixed, under some 
circumstances it can evolve towards a pure state again. 

If instead of Ucnot(0) the interaction were Uswap(@), the noise would affect the system 
in an different way. In this case the system will be left in one of the states 



Po = — 
Po 

1 

Pi = — 
Pi 



wo\i>o)(i>o\ +wi\a\ 2 sin 2 6\l)(l\ 
Wl |^i)^i|+^o|/?| 2 sin 2 ^|0)(0| 



(112) 



where 



and 



|^ ) = a \0) +(3cos6\l) , 

= acos#|0) +(3e- 2i6 \l) , 

Po = Wq cos 2 6 + | or | 2 sin 2 9 — 1 — p\ 



113) 



(114) 



In this case, rather than evolving eventually towards one of the pure states |0), |1), the system 
will evolve towards the state u>o|0)(0| +wi|l)(l|, i.e., it will become identical with the mixed 
state of the environment bits. This is quite similar to the process of thermalization, by 
which a system interacting with a heat bath evolves towards equilibrium. 



C. Trajectories with generalized measurements 

Now we will assume that both the system and environment q-bits begin in a pure state, 
but that the measurement performed on the environment bit is no longer a projective mea- 
surement, but instead a POVM which provides only partial information about the state of 
the environment bit. 

Let us consider a system where the system and environment q-bits interact via 
^s(^)^cnot(^), the environment bits start initially in the state |0), and the system is in 
state a|0) + /3|1). After the system and environment interact, the environment is measured 
with a POVM using operators 

E = g|0)(0| + (l-g)|l)(l|=i 2 , 
^ 1 = (l-g)|0)(0|+g|l)(l|=i 2 , 

1 = v /^|0>(0| + v /T^|l)(l|, 

1 1 = v /T^|o)(o|+Vg|i)(i| ■ (H5) 
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Assume that 1/2 < q < 1. A measurement result then favors the environment q-bit being 
in state |0), and a result 1 favors |1), but neither result is conclusive. 

After the system and environment interact, the probabilities of the two outcomes are 

Po = q(\a\ 2 + \f3\ 2 cos 2 6) + (I - q)\{3\ 2 sin 2 6 = I - Pl . (116) 

If we assume S < 1, tracing out the environment leaves the system in one of the mixed 

states 

p « (l - (1 - q)\(3\ 2 9 2 /q) mm + ((1 - q)\(3\ 2 9 2 /q) |1)(1| , 



Pi 



(l - q\(3\ 2 6 2 /(l - q)) IVWI + (q\(3\ 2 e 2 /(l - qj) |1)(1| , (117) 



with the updated state 

= a (l + \P\ 2 9 2 /2)\0) + /3(1 - \a\ 2 9 2 /2)\l) . (118) 

That is, the system state remains entangled with the environment q-bit. 

This environmental interaction and measurement leaves the system in a state of the form 

p 3 = w\if/)(il/\ + Q.-w)\l)(l\ . (119) 

If a state in that form interacts with another environment bit, on which the same POVM 
is performed, the system state will remain in the same form. After n interactions, the state 
will have the form ( |111|) , with \%p n ) — ► |0) as n —>■ oo, and the state tending towards either 
|0)(0| or |1)(1| at long times. Thus, the system evolution looks quite similar to the case when 
the environment is initially in a mixed state and a projective measurement is performed. 
However, in the case of a mixed environment the system has no entanglement with the 
environment, and is in a mixed state due to noise; in the case of incomplete measurements, 
the system is entangled with the environment, and the system and environment as a whole 
remain in a pure state. 

In this case as well, one can describe the evolution of the system by a stochastic master 
equation. This has exactly the same form as ( 110f) , with u>o,i replaced by q and 1 — q. Unlike 



the previous case, however, the system state is not mixed because of classical noise entering 
from the environment, but rather because it remains entangled with the environment. 

The POVM just considered represents a measurement with imperfect discrimination. 
Let's consider a different way that a measurement may be incomplete. Consider the following 
positive operators: 



E 


= g|0><0| = 


A 2 


E x 


= <?|1><1| = 


Al, 


E 2 


= (1-9)1 


= Aj 


A 


= V#><o| 


•) 


h 


= ^ll>(ll 






= \/i-?i 





(120) 
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These operators Ei sum up to the identity, and for < q < 1 form a POVM. Eq and 
Ei are proportional to orthogonal projectors; E2, by contrast, leaves the state completely 
unchanged. This is like a measurement made with a detector of efficiency q; there is a 
probability 1 — q that the detector will fail to register anything at all, and hence gives no 
information about the state. 

Suppose the system q-bit is originally in a state of form 

p s = w \^)^\+w 1 \l}(l\ . (121) 

(where an initial pure state would have w = 1). Then a measurement result of has 
probability 

p = q [w (\a\ 2 + \(3\ 2 cos 2 9) + w 1 cos 2 e\ « q[l - (1 - w \a\ 2 )9 2 ] , (122) 
and leaves the system in a mixed state of the same form 

po « w (l + Wl \a\ 2 ^)W)W\ + M 1 - w \a\ 2 9 2 )\l)(l\ ; (123) 



the result 1 has probability 



V\ = q 



wo|/3| 2 sin 2 + wisin 2 « q{\ - w \a\ 2 )6 2 , (124) 



and leaves the system in the state p\ = |1)(1|; and the result (or rather, nonresult) 2 has 
probability P2 = 1 — q and leaves the system in the mixed state 

P2 » ™ (i - \P\ 2 o 2 )W)W\ + K + w \p\ 2 e 2 )\\){\\ . (125) 

If the evolution began in a pure state, it will remain in a pure state unless outcome 2 occurs; 
but once this happens, the state will tend to remain mixed until either a result 1 occurs 
(in which case the system will be in state |1) thenceforth) or, in the limit of large n, the 
system approaches the state |0). The measurement process works much as before, but the 
null results slow the convergence; if the efficiency q is low, it can take far longer for the state 
to settle down to either |0) or |1). 



X. QUANTUM TRAJECTORIES, CONSISTENT 
HISTORIES, AND COLLAPSE MODELS 

As is clear from the previous sections, the formalism of quantum trajectories calls on 
nothing more than standard quantum mechanics, and as framed above is in no way an alter- 
native theory or interpretation. Everything can be described solely in terms of measurements 
and unitary transformations, the building blocks of the usual Copenhagen interpretation. 

However, many people have expressed dissatisfaction with the standard interpretation 
over the years, usually due to the role of measurement as a fundamental building block of the 
theory. Measuring devices are large, complicated things, very far from elementary objects; 
what exactly constitutes a measurement is never defined; and the use of classical mechanics 
to describe the states of measurement devices is not justified. Presumably the individual 
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atoms, electrons, photons, etc., which make up a detector can themselves be described by 
quantum mechanics. If this is carried to its logical conclusion, however, and a Schrodinger 
equation is constructed for the measurement process, one obtains not classical behavior, but 
rather giant macroscopic superpositions such as the famous Schrodinger's cat paradoxes. 

Three major approaches have been followed in tackling this problem. In the first, "hidden 
variables" are introduced which determine the outcome of measurements, so that quantum 
mechanics is a partial description of an underlying deterministic theory. Unfortunately, as 
John Bell showed@, such a hidden variables theory (such as that of de Broglie and Bohmil) 
must be nonlocal to give the same predictions as quantum mechanics. 

The second approach is to modify quantum mechanics to get rid of the unwanted macro- 
scopic superpositions, while retaining the usual quantum results on small scales. Theories of 
this nature have been proposed by Pearle0, by Ghirardi, Rimini and Weber0, by Diosill, by 
Gisini and Percival0, and by Penrose^, among others; these theories are commonly called 
"collapse models." 

These alternative theories usually replace the Schrodinger equation with a new, stochastic 
equation, which includes a mechanism to produce wave function collapse onto some preferred 
basis. Such stochastic Schrodinger equations often strongly resemble those produced by a 
suitable quantum trajectory description, in which the system interacts continuously with 
an environment which is repeatedly measured. While I do not know if all such alternative 
theories are equivalent to a quantum trajectory description, certainly a large class of them 
is. 

Because these theories are not equivalent to standard quantum mechanics, in principle 
they can be distinguished experimentally. Unfortunately, most such theories require experi- 
mental sensitivities well beyond current technology, though Percivaffi has made interesting; 
proposals of possible near-term tests. 

Theoretically there are also problems in producing relativistic versions of such theories. 
Most such models produce collapse onto localized variables, such as position states. Such 
localization, however, does not in general commute with the Hamiltonian, and so such the- 
ories almost always violate energy coaservation-sometimes spectacularly soS Another 
class of models produces a collapse onto states of definite energyillllll, rather than local- 
ized variables. It is not clear, however, that such collapse models eliminate macroscopic 
superpositions in either the short or long term. Certainly the end result, a universe in an 
energy eigenstate — a static, unchanging state with no dynamics — does not greatly resemble 
the world that we perceive, though in the context of an expanding universe even this is 
unclear. 

The last approach is to retain the usual quantum theory, but to eliminate measurement 
as a fundamental concept, finding some other interpretation for the predicted probabilities. 
While there are many interpretations that follow this approach, the one that is most closely 
tied to quantum trajectories is the decoherent (or consistent) histories formalism of Griffiths, 
Omnes, Gell-Mann and Hartleiffl. In this formalism, probabilities are assigned to histories 
of events rather than measurement outcomes at a single time. These can be grouped into 
sets of mutually exclusive, exhaustive histories whose probabilities sum to 1. However, not 
all histories can be assigned probabilities under this interpretation; only histories which lie 
in sets which are consistent, that is, whose histories do not exhibit interference with each 
other, and hence obey the usual classical probability sum rules. 



31 



Each set is basically a choice of description for the quantum system. For the models 
considered in this paper, the quantum trajectories correspond to histories in such a con- 
sistent set. The probabilities of the histories in the set exactly equal the probabilities of 
the measurement outcomes corresponding to a given trajectory. This equivalence has been 
shown between quantum trajectories and consistent sets for certain more realistic systems, 
aswel&H. 

For a given quantum system, there can be multiple consistent descriptions which are 
incompatible with each other; that is, unlike in classical physics, these descriptions cannot be 
combined into a single, more finely-grained description. In quantum trajectories, different 
unravelings of the same evolution correspond to such incompatible descriptions. In both 
cases, this is an example of the complementarity of quantum mechanics. 

We see, then, that while quantum trajectories can be straightforwardly defined in terms 
of standard quantum theory when the environment is subjected to repeated measurements, 
even in the absence of such measurements there is an interpretation of the trajectories in 
terms of decoherent histories. Because the consistency conditions guarantee that the prob- 
ability sum rules are obeyed, one can therefore use quantum trajectories as a calculational 
tool even in cases where no actual measurements take place. 



XI. CONCLUSIONS 

In this paper I have presented a simple model of a system and environment consisting 
solely of quantum bits, using no more than single-bit measurements and one- and two-bit 
unitary transformations. The simplicity of this model makes it particularly suitable for 
demonstrating the properties of quantum trajectories. We see that these require no more 
than the usual quantum formalism, though they can, of course, also be applied in cases going 
beyond standard quantum mechanics, such as collapse models or the decoherent histories 
interpretation. 

Quantum trajectories can often simplify the description of an open quantum system in 
terms of a stochastically evolving pure state rather than a density matrix. While for the 
q-bit models of this paper there is no great advantage in doing so, for more complicated 
systems this can often make a tremendous difference!. 

In this model it is also possible to straightforwardly quantify the flow of information 
between system and environment, in the form of entanglement and noise, as well as the 
information acquired in measurements and the decrease in entropy of the system. For 
realistic systems and environments this may be more difficult to do analytically, but one can 
still form a qualitative picture of the information flow in decoherence and measurement. 

I hope that this model will clear up much of the confusion that currently surrounds the 
theory of quantum trajectories. Quantum trajectories are a very useful formal and numerical 
tool, and should find applications well outside their current niches in quantum optics and 
measurement theory. 
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APPENDIX: SOLUTIONS TO EXERCISES 

Exercise 1. Any unitary operator can be written U = exp(iH) where H is some 
Hermitian operator. Any operator on a single q-bit can be written in terms of the basis 

H = c l + Ci& x + c 2 <7 y + c 3 a z . 

For H = W to be Hermitian, the q must all be real. Choose a convenient parametrization 
ci,2,3 = <f>n x , y ,z, where n 2 + n 2 + n 2 — 1. The expression for U is then 

U = exp(icol + i(pn ■ a) = exp(icol) exp(i(pn ■ a) , 

where the Co term factors out because 1 commutes with everything else. This factor just 
contributes a global phase change, and hence can be neglected. The second exponential can 
be expanded in a power series to yield 

oo 

exp(i(ftn ■ a) = ^ (i<fm ■ a) rn jm\ . 

m=0 

Note now that 

(n . a) 2 = (nl + n 2 y + n 2 z )l + n x n y (a x a y + a y a x ) 

+n y n z (a y a z + o z a y ) + n z n x (a z cr x + o x a z ) 
= i 

where we've used the facts that a x ,a y , a z anticommute and that a 2 — a 2 — a 2 — 1. Substi- 
tuting this into the power series yields 

oo oo 

U=J2 (i<P) 2m (n ■ B) 2m /(2m)\ + (ft ■ &) ^ (^) 2m+1 (n • B) 2m /(2m + 1)! 

m=0 m=0 

oo oo 

= i W>) 2 "7(2m)! + (n ■ 3) £ (#) 2m+1 /(2m + 1)! 

m=0 m=0 

= i cos + i(n ■ a) sin , 
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which was to be shown. 

Exercise 2. The proof of this is just the same as in Exercise 1. Simply expand the 
exponential in powers of its argument 



U(Q) = exp(-zW) = £ (-i6U) m / 



ml . 



m=0 



Because U 2 = WU = 1, all the odd powers will be proportional to U and the even terms to 
1. Collecting the even and odd powers separately gives the result 

U(6) = Ucos6 ~ ill sin6 , 

just as in Exercise 1 above. 

Exercise 3. The average information gain from a measurement is given by the Shannon 
entropy (|26|). 



St, 



-po log 2 Po - Pi log 2 Pi 



In this case, the relevant probabilities are po = 1 — e|/9| 2 and p\ = e\/3\ 2 . If we rewrite 
log 2 x = lax/ In 2, we can use the expansion ln(l + 5) = 5 — 5 2 /2 + <5 3 /3 — • • • to get 



Smeas~e|/3| 2 (l/ln2-log 2 e| 



Exercise 4. Because p = is an Hermitian operator, we can always find an orthonormal 
basis which diagonalizes it: 





f\l 





■ 


' \ 







A 2 


■ 


• 


p = 








A 3 ■ 


■ 




V o 





• 


• aJ 



where the {Aj} are all real. The positivity of p requires that Aj > 0, and Trp = 1 implies 
Y, Aj = 1. Putting all this together, we find that the von Neumann entropy becomes 



-Tr{p log 2 p} = - J2 Ai log 2 A; 



8=1 



The right-hand side of this equation is exactly the expression for a Shannon entropy, as in 



(2"Sp. For an n-outcome measurement, the Shannon entropy has a single global maximum 
at Ai = A 2 = ■ ■ ■ = A n = 1/n. If we substitute these values into the matrix expression for p 
above, we find that p = 1/n = 1/Trl. 
Exercise 5. 

Y,dldk = Y, x kPkjPl f A],A 3 

k jj'k 

jj' 
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33' 

^Tr cnv {{A),A j ® B^Bj) (i s ® |£>(£|)} 



= Tr env {£t£(ls® 

= IV env {i s ® \E)(E\} 

= U- 

Exercise 6. A very direct way of showing this result is to use the simple approximation 
1 - 9 2 ps exp(-# 2 ) for 9 < 1. This gives: 



n na-i^i^ 

n=0 
oo 

n 

n=0 
oo 

n 



n=0 n=0 

00 |a| 2 + (l-# 2 )|/3| 2 exp(-nfl 2 ) 



a\ 2 + \(3\ 2 exp(-n6 2 ) 



00 'a| 2 + |/3| 2 exp(-(n + 1)# 2 



n=0 



a| 2 + |/3| 2 exp(-n# 2 ) 



lim 

71— >0O 



\a\ 


2 + \(3 


| 2 exp(-(n + l)# 2 ) 




a\ 


2 + \f3\ 2 



2 



A slightly more indirect approach is to take the natural logarithm of the first line of the 
above equation, and use the approximation ln(l + x) ~ x for i < 1. This gives 



OO 



In J] Pd(n) w In n (1 ~ WO 2 ) « - E IA»| 2 * 2 ■ 

n=0 n=0 n=0 

By approximating the sum as an integral and then exponentiating both sides, one straight- 
forwardly gets the same result. 

Exercise 7. Using the definition fl87p, we see that 



E^p4 = ^(|0)(0|+e^|l)(l|)p(|0)(0|+e^|l)(l|) 
+ (|0)(0|+e^|l)(l|)p(|0)(0|+e^|l)(l|) 



2 



By rearranging and combining terms, we get 



Ei fe p4 = ^((|O)(O| + (cos0-ism0)il)(l|> 
x(|0)(0| + (cos0 + isin0)|l)(l|) 
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+(|0)(0| + (cos0 + zsin0) 1 |)p 

x(|0)(0| + (cos6 -ism6)\l)(l\) 

= (|0)(0| + cos^|l)(l|)p(|0)(0| + co^|l)(l|) 
+ sin 2 0|l>(l|p|l>(l| , 

where the last equality is just A pA + AipA| using the definition (|81|). 

Exercise 8. The procedure here is just the same as in Exercise 7. We write the new 
state in terms of the operators ( |96l) and rearrange the terms to get the operators (|Sll). 

k 

= ^((|0)(0| + (cos0 - sin0)|l)(l|)p(|O)(O| + (cos6 - sin (9) 1 1> <1|) 

+ (|0)(0| + (cos0 + sin0)|l)(l|)p(|O)(O| + (cos0 + sin 6>) 1 1> ( 1 1 ) 

= |0)<0|p|0}<0| + 2cos0(|O>(O|p|l><l| + |l)(l|p|0)(0|) 

+ cos 2 0|l>(l|p|l>(l| + sin 2 0|l)(l|p|l)(l| 
= (|O)(O| + cos0|l)(l|)p(|O)(O|+cos^|l)(l|) 

+ (sin0|l)(l|)p(sin0|l)(l|) . 

Of course, the evolutions given by (|8l|) and (|7]) have already been shown to be equivalent, 
so the equivalence of (|96|) and (^7]) follows at once. 

Exercise 9. In this measurement, the system goes from the maximally mixed state 
p = 1/2 to one of the two states p ,i given by ( |103| ). The von Neumann entropy of the 
maximally mixed state is 

fl ll 1, 1 1, 1 

- Tr \2 l0g2 2/ = -2 l0g2 2-2 l0g2 2 = 1 - 

Because the states po,i are diagonal, it is straightforward to evaluate their von Neumann 
entropy as well; we find it to be 

-Tr{p ,ilog 2 po,i} = -p log 2 Po -pilog 2 pi = S meaa , 

where po,i = (1 ± 6>)/2. So we see that 

Smeas = 1 ~ ~ log 2 (l + 9) - ^ l g 2 (l - 0) « 1 - 2 /(2 hi 2) , 

which implies that AS* ~ 2 /(21n2) and AS/S mcas — > as 9 — > 0, as required. The entropy 
^meas of the measurement is nearly a full bit, but only a tiny part of it represents information 
about the system. 
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FIGURES 




Environment 



FIG. 1. This is a schematic diagram of the type of model used in this paper. The system is a 
single two- level system, or q-bit. It interacts briefly with a series of environment q-bits which pass 
by with average time intervals of St, and the environment bits are subsequently measured. 
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Cavity 



Incoming radiation modes (vacuum) 




Fully-silvered Partially-silvered Outgoing radiation modes (excited 



FIG. 2. Here is a particular physical system of which the model of this paper is an idealization. 
The system is an electromagnetic field mode in a cavity with a partially-silvered mirror. There 
is on average less than a single photon in the mode at a time, so only the lowest two occupation 
levels contribute significantly to the dynamics. The external electromagnetic field is described in 
terms of incoming and outgoing wave packets of width 5t, each of which reflects on the outside of 
the partially-silvered mirror. The incoming packets are all in the vacuum state, but upon reflection 
they may absorb a photon emitted from the cavity. 
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(B) 

FIG. 3. These two figures show a variety of quantum jump trajectories for different initial 
states, plotting the squared amplitude |/3| 2 for the system to be in state |1) versus the time. Figure 
a) assumes an interaction £7cnot with the environment, while b) assumes [/swap- In both cases, 
\(3\ 2 decays towards zero. Some trajectories just decay smoothly to 0, while others undergo a jump. 
In figure a) this jump is to 1, while in b) they jump to 0, as in spontaneous emission. Note that 
the decays of |/3| 2 is not a simple exponential, due to the nonlinearity of the trajectory equations. 
As \(3\ 2 — ► 0, however, it approaches an exponential. 
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FIG. 4. This figure plots |/3| 2 against time for the diffusive equation (p5|), with several different 
initial states and realizations of the stochastic term. We see that when the trajectories approach 
or 1 they tend to remain there, while in between they diffuse freely. 
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